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Abstract. We derive the exact Casimir-Polder potential for a polarizable microparticle inside an ideal 
metal cylindrical shell using the Green function method. The exact Casimir-Polder potential for a particle 
outside a shell, obtained recently by using the Hamiltonian approach, is rederived and confirmed. The exact 
quantum field theoretical result is compared with that obtained using the proximity force approximation 
and a very good agreement is demonstrated at separations below 0.1/?, where R is the radius of the cylinder. 
The developed methods are applicable in the theory of topological defects. 



1 Introduction 

The Casimir-Polder and Casimir forces act between a po- 
larizable microparticle and a surface and between two sur- 
faces, respectively [U[l]. Both forces are of quantmn na- 
ture and originate from zero-point oscillations of quantized 
fields. The Casimir force due to nontrivial spatial topology 
plays an important role in the compactification of extra 
dimensions in Kaluza-Klein theories [3] , in cosmology 4^ , 
and in the theory of topological defects [5] . With advances 
in nanotechnology and successes in the production of ul- 
tracold atoms it became possible to measure the Casimir 
and Casimir-Polder force with increased precision (see [6l 
[7] for a review). Interestingly, some of the fundamental re- 
sults on the Casimir effect have already found topical ap- 
plications. As an example, the Casimir energy for isolated 
cylindrical and toroidal carbon nanotubes was calculated 
in OE] basing on the theory of topological defects. 

Many calculations of the Casimir and other interac- 
tions between curved surfaces were performed by means of 
the so-called proximity force approximation (PFA) [1211131 
1141115) . This is an approximate method which replaces in- 
finitesimal elements of curved surfaces with pairs of plane 
parallel plates and uses advantages of a known exact solu- 
tion for the plane parallel geometry. The PFA is commonly 
used for comparison between experimental results and the- 
ory [16]. This method is subject to some restrictions and 
possesses only a limited accuracy. Because of this, a more 
fundamental quantum field theoretical approach is highly 
desirable. Such an approach was developed, for instance, 
for an ideal metal cylinder above an ideal metal plane 
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|10) . For this configuration, the asymptotic expansion of 
the exact result in the limit of short separations was found 
and compared with the PFA [TT]. As one more configura- 
tion allowing the fundamental description based on first 
principles of quantum field theory one could consider the 
interaction potential of a microparticle with an ideal metal 
cylindrical shell. This configuration has long been studied 
(see related references in 17 ). Recently, the exact interac- 
tion potential for a microparticle external to a cylindrical 
shell has been found in |17| using the Hamiltonian ap- 
proach. It was also suggested [TH] to use Rydberg atoms 
out of thermal equilibrium inside a cylindrical cavity for 
the observation of resonant Casimir-Polder interactions. 

Here, we derive the exact expression for the Casimir- 
Polder potential between a polarizable microparticle ei- 
ther external or internal to an ideal metal cylindrical shell 
using the Green tensor of the electromagnetic field in a 
cylindrical configuration. For a microparticle external to 
a cylindrical shell the obtained results are shown to be 
in agreement with those derived in |17j with the help of 
another method. We compare the exact particle-cylinder 
potential with the approximate result obtained using the 
PFA. This allows one to determine rigorously the appli- 
cation region of the PFA in one more of rare cases in the 
Casimir physics where the exact quantum field theoretical 
solution is available. The role of the dynamic polarizability 
of a particle is investigated. The suggested exact approach 
based on the use of the Green tensor allows generalization 
for different boundary conditions. 

The structure of this paper is as follows. In Sect. 2 
we derive the Green tensor for an exterior region of a 
cylindrical shell. In Sect. 3 the same is done for an inte- 
rior region. Section 4 is devoted to the derivation of the 
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exact Casimir-Polder potential for a particle situated ei- 
ther in exterior or interior regions of the cylindrical shell. 
In Sect. 5 the role of dynamic effects is considered and 
the comparison between the exact and the PFA results is 
performed. Section 6 contains computational results for a 
particle inside a cylindrical shell. In Sect. 7 the reader will 
find our conclusions and discussions. 



2 Green tensor for an exterior of ideal metal 
cylindrical shell 

We consider the Casimir-Polder interaction between a mi- 
croparticle (an atom or a molecule) and an ideal metal 
cylinder. The corresponding interaction energy is given 
by the expression (see, for instance, [TO] ) 



(1) 



where aji(i£,) is the polarizability of a particle situated at 
a point r calculated along the imaginary frequency axis, 
J, Z = 1, 2, 3, and 
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Here, x = {t,r), x' = (t',r'), r = i - i'. In ©, Gj;^(a;,a;') 
is the part of the retarded Green tensor Gji{x, x') induced 
by the cylindrical shell: 



G,i{x,x')^Gf{x,x') + G'^;'^{x,x'), 



(3) 



where G^^^{x,x') is the Green tensor in free space. The 
retarded Green tensor is defined as 

hGji{x,x') = -ie{T){E,{x)Ei{x') - Ei{x')Ej{x)), (4) 

where 9{x) is the unit-step function, Ei{x) is the operator 
of the i component of an electric field, and the angular 
brackets mean the vacuum expectation value. 

For the evaluation of the vacuum expectation value in 
([U we use the mode sum formula 

Gfl{x,x') = -i6{t)Y}E„A^)Ki{^') - EM)El^{x)], 

a. 

(5) 

where {£0(0;), E* (x)} is the complete set of eigenfunc- 
tions for the electric field in the region outside a con- 
ducting cylindrical shell, specified by the collective index 
a — (7,m, fc) (see below), and the asterisk stands for a 
complex conjugate. For the geometry under consideration 
we have two different classes of eigenfunctions correspond- 
ing to the cylindrical waves of the transverse magnetic 
(TM) and transverse electric (TE) types. In cylindrical co- 
ordinates (r, 0, z), the eigenfunctions, satisfying the ideal 
metal boundary condition on the surface of a cylinder of 
radius i?, have the form 



(6) 



where m = 0,±1,±2,..., —00 < k < 00, and A = 1,2 
correspond to the TM and TE waves, respectively. For 
separate components of the electric field the radial func- 
tions in ^ have the form 
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where here and below oj = cy/^-f+k^, < 7 < cxd, the 
prime means the derivative with respect to the argument, 
and the following notations are introduced: 

gi,a{ir) = J\,n\{ir)Y\,^l{'^R) - ^Iml (7'')'^|m| (7^), 

92,ahr) = Ji,n\{ir)Y;^\{lR) - >|™|(704n|(T^)- (8) 

In equation J,y{x) and Y^{x) are the Bessel and Neu- 
mann functions. The normalization coefhcient in (|6]) is 
given by the expression 



(9) 



Substituting the eigenfunctions ([B]) into the mode-sum 
([S]), we find for the Green tensor 



G,,(x,x') = -*^(r)^ ^ Pi, 



(10) 



a A=l,2 



To.Ei']{r)Elf{r') - T:£;2(r')i?S*M' 



Here, the summation over the collective index a means 
the summation over m and integrations over k and 7, and 

Ta — exp{imA(j) + ikAz — ilot), 
A(j)^(j)' - <j), Az^ z' ~ z. 

In order to extract from ([TU)) the cylinder- induced part, 
we subtract from the right-hand side of this equation the 
respective expression for the bulk without a cylindrical 
shell. The latter is given by (ITUl) , where should be 

replaced with c^/{2it^uj) and the respective field compo- 
nents are given by ([7]) where the functions g\^a{l'r) are 
replaced with J\m\(^r). As a result, the cylinder-induced 
part of the Green tensor can be presented in the from 



G 



70; 



where the following notation is used: 
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Here, the quantities E^j'^\r) are given by ([7]), where the 
functions gx.ai'^r) are replaced by the Hankel functions 
H^^^i^r). Tilde in ^ and ^ stands for the complex 
conjugation which is, however, not applied on the func- 



tions _fff*l(7r) and H}''''.{jR). For example 



By taking into account ([TT|) . for the spectral compo- 
nent of the Green tensor in ((1} one finds 



a ' s=l,2 I 



E 



E' 



a,l 
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uj — LJ + 



(13) 
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a; — 



uj + 



Using the expressions for E'^'J^\r), for the off-diagonal 
components we have 

Gj)(r,r;zO = -GgHr,r;zO. 

Since the tensor aji{i^) is symmetric, we conclude that the 
off-diagonal components do not contribute to the Casimir- 
Polder interaction energy. 

For the diagonal components, in ((13)) we rotate the 
integration contour in the complex 7-plane by the angle 
tt/2 for the s = 1 term and by the angle — 7r/2 for the s — 2 
term. The integrand has poles at 7 = ±i^yk^+^^/(? and 

we assume that they are avoided by small semicircles Cp^^ 
of radius p in the right half-plane. As a result the integrals 
over 7 for the s = 1 and s = 2 terms are presented as the 
sum of integrals over the interval (0, ±iy^k^ -I- ^^/c^ =F«p), 

over C^=^' and over the interval (zti-^fc^ -|- £^'^/c'^ ± jp, 00). 
It can be seen that the integrals over the two intervals are 
cancelled for the s = 1 and s = 2 terms. Evaluating the 

integrals along C^*"* , for the diagonal components of the 
Green tensor we find the following expressions 



c2 K'„^{I3R) /32r2 



Km{f3R) 
Gg(r,r;*0 = -E / 

m=0 "^^ 

-k 



00 r (2 Ti 



c2 K'^m) 



2 I^m) m^KliPr) 



4 °° / /■= 
Gg(r,r;ze) = — ^ / 



(14) 



where Im{z) and Km{z) are the modified Bessel functions 
and we introduced the notation /3 = y/k'^ + jc?-. The 
prime near the summation sign means that the term with 
m = is multiplied by 1/2. 



3 Green tensor for an interior of ideal metal 
cylindrical shell 

In this section we consider the region inside a cylindrical 
shell. The corresponding eigenfunctions have the form ([5]), 

where the functions Eq^^(7') are given by expressions ^ 
with the replacement gx^ai'yf) — ^ J|„j|(7r) and for the 
normalization coefficient one has 



ttRoj'j 



jmx) + {l-iyyx^)J^{x)' 
(15) 

From the boundary conditions on the cylindrical bound- 
ary with radius R, we can see that the eigenvalues for the 
quantum number 7 are the roots of the equation J|m| (7-R) = 
for TM waves and the roots of the equation J'^^_^^ (jR) — 
for TE waves. The corresponding eigenvalues ^R we de- 
note jR = jm]n, n = 1,2,..., where, as before, A = 1, 2 
correspond to the TM and TE waves. The mode sum for 
the Green tensor is given by (jlOp , where now the summa- 
tion over the collective index a means the summation over 
m and n, and the integration is performed over k. For the 
spectral component of the Green tensor one finds 



G,,(r,r;zO = --E E 



2] ml {jm}n) 



.(A) (A) 

A— 1 2 J'm,n^m,n 



(16) 



Ei,>ir)Elfir) ^ E^'J (r) E^^j* (r) 



(A) 



where UJm^n = C\J jm^n/R"^ + fc^- 

For the summation of the series over n in (|16p we apply 
the formula pU] 

E^i™iOtl)/(ji'l) = o / d^fi^) (17) 

n=l •'0 

\m\ V / 



where 



and 



(^) = {z), Ji^i (z) = J|',„| (z), 



yEi^««I.=.../(^)- ,(a). . 



(2), 



Formula (jl7p is valid for functions f{z) analytic in the 
right half-plane Rez > with poles ±ii/p, Up > 0, p — 
1,2,..., and satisfying the condition 



/(ye 



TTi/2\ 



„2|m|Tri J(ye~'^*/^ 



4 



V. B. Bezerra et al.: The exact Casimir-Polder potential and the proximity force approximation 



and the condition 



\fix + iy)\<e{x)e''y, h<2, 

for ?/ — cxD, with e{x) — )■ for x ^ oo. For the function 
f{z), corresponding to the components of the Green tensor 
(fT6|) . the poles on the imaginary axis correspond to yp = 



+ ^2 and the above conditions are satisfied \fr < R. 
For example, in the contribution of the TM waves (A = 1) 
to the rr-component of the Green tensor the function f{z) 
has the form 

f{z) = zJ{l^{zr/R)/[z^ + R^{e+e)]- 

The part of the Green tensor corresponding to the first 
term on the right-hand side of P7)) is the Green tensor for 
the free space. The boundary-induced part of the Green 
tensor corresponds to the second term. It can be seen that 
the diagonal components of the boundary-induced part in 
a cylindrical cavity are obtained from (|Ti)) by the replace- 
ments Ira ^ Km- 



4 The Casimir-Polder potential 

We first consider a particle outside of a cylindrical shell. 
Substituting ([TH) into ((H we obtain the exact formula for 
the Casimir-Polder interaction energy in the general case 
of an anisotropic polarizability 



U{v) 



2h 



d(3 



c2 Ki^im > 



2 Im{l3R) 2 .O N I 



(18) 



where am ct<f,<f,, and azz are the diagonal components of 
an atomic (molecular) polarizability. 

In the isotropic case the exact Casimir-Polder energy 
(fTS]) simplifies to 



U{r) 



2h 



m— 



^^"(*^)y^//^y^2_^2/c2 



G™(/3r) 



with the notation G„i{z) — K'^iz) + {m/zf' K^iz). 



(19) 



It can be easily seen that (fT5)) coincides with expres- 
sion (17) [with notations (27)-(29)] of [T^, where the ex- 
act Casimir-Polder potential for an atom outside an ideal 
metal cylinder was obtained. To see this we use the fol- 
lowing expressions for the components of atomic polariz- 
ability 

Here, \i) is the state of an atom with the energy Ei, Eji = 
Ej—Ei, and ijl is the atomic dipole moment operator. Now 
we substitute ((20)) into (fT8|) and return to the integration 
variable k: 



^ J rn=0 Jo 



dk 



1 



El + h-^e 



c2 K^{I3R) /32r2 

e I'ram) 



c2 K'^{PR) 



KiWr) - k^ 



2 Imi^R) r^>2rg . 



(21) 



Next we introduce polar coordinates in the plane (^/c, k): 
^/c — P cos 9, k = /3 sin 0. This allows us to rearrange (PT|) 
into the form 



i m=0 •^'J •^'J J 



P^d9 



+ h^c^/S^ cos2 



cos 



2 



(22) 



cos^e4^<(/3r) 



lUPR) m^KliPr) 



1 2 ImiPR) 

' KmiPR) 



Evaluating the integrals over 9, we obtain from the 
following result: 




d/3/3 



(23) 



+ \lJ.r\ 



E, 



I El + h^c^p^ 



1 \E 



urn) m^KUM 



El + K^c^p^ - Ej, 



I,n{PR) 
KmiPR) 
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The latter coincides with the result given by equations 
(17), (27)-(29) of tli (note that in [T7] the International 
system of units is used) . 

For an anisotropic particle inside a cylindrical shell the 
Casimir-Polder potential is obtained by the substitution 
of the Green tensor (jl6p in Q , which results in 



dl3 



X < a. 



c2 /?V2 

•2 ^ 



(24) 



The Casimir-Polder energy for an isotropic particle inside 
a cylindrical shell is given by the expression 



U{r) 



2h 



m— 



h/c v/3 - r/c-^ 



\e K'^m) p KmWR) 



(25) 



where Fjn{z) — Im{z) + {m/z)^ ^mi^)- We emphasize that 
the results (j24l) and (j25p cannot be obtained from the for- 
malism [21 developed for an atom inside a hollow dielec- 
tric cylinder of finite permittivity. 



5 The role of dynamic effects and the 
comparison with the proximity force 
approximation 

The following analysis is given for the expression (IT9|) . 
i.e., for an isotropic polarizable particle spaced outside a 
cylindrical shell. 

Let us denote by a the separation distance between a 
particle and a surface of the cylinder, so that a = r — R. 
At small distances, a/R <^ 1, the dominant contribution 
to Eq. (jl9p comes from large values of m and we can use 
the uniform asymptotic expansions for the modified Bessel 
functions for large values of m. After calculations, in two 



leading orders in a/_R, we find 
U{a) 



167ra4 



dCa(jWcC) {(2 + 2C + C^)e-f (26) 

[(2 + 2C + C')e-^ + 2C^r(o, C) - C'/^(-2, C)] } , 

where uJc = c/ (2a), C = ■C/'^c, and r{b, C,) is the incomplete 
gamma function. The term of zeroth order in a/i? in (|26p 
coincides with the Casimir-Polder potential for an atom 
near an ideal metal plate (see equation (16.27) in |22)V 

The oscillator model for the dynamic polarizability of 
a particle 



9] 



(27) 



where gj are the oscillator strengths and ujj are the eigen- 
frequencies, works well over a wide range of separations 
[231124] . Substituting ^ into (HH) and evaluating the in- 
tegrals over ^ explicitly, one finds 



j 3 m— 

G,„(/3r) /™(/3i?) 



I'miPR) 



'^^^\K'^m) 



(28) 



,(/3)[s,(/3) + l] KraiPR) 



,(/3) s0) + l 



(p-pp- jiij^- + 1. The main contribution to 



where Sj(/3) e 

the /3-integral in ([25)) comes from /3 < 1/a. If a c/ojj, 
the ratio c/3/wj is small and from ([25]) we obtain 



OO r, 

t/(-)(o = ^c^E7 

^-n "'0 



ImWR) 



[G„,(/3r)+2i^2,(/3r)] 



(29) 



where uq — gj /uj"^ is the static polarizability of a par- 
ticle. In the static limit, assuming a/R -C 1, from (l26l) we 
find 

f/(--)(a)«[/o(a)(l-|^), (30) 

where C/o(fl) = — 3?icao/(87ra'*) is the famous Casimir- 
Polder result for the interaction energy between a particle 
and an ideal metal plane [T]. 

In Fig. [T] we plot the exact Casimir-Polder energy (l28t 
for a particle outside a shell normalized for Uq as a func- 
tion oi a/R using the single-oscillator model [ ([?7|) with 
wi = Wyi]. Different hues are plotted for R/Xa — 1, 2.5, 
and 5 (lines 1, 2, and 3 respectively), where Xa — 2Trc/u!A- 
The dashed line corresponds to the static limit ([2U)) . For 
instance, for an atom of metastable He* one has [53] uja = 
1.18 eV and the cylinder radius corresponding to R/Xa = 
1 is equal to 0.95 /j,m. As can be seen in Fig. [I] with in- 
creasing R/Xa the Casimir-Polder energy approaches to 
the static limit. Note that although lines 1,2, and 3 around 
a/R — demonstrate the same distance dependence (the 
inverse third power of separation) as does the nonrela- 
tivistic limit (i.e., the van der Waals force), they are quite 
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1.1 - 
1.0- 



0.9- \ 




0.4l 



0.5 1.0 1.5 2.0 

a/R 

Fig. 1. The normaUzed exact Casimir-Polder potential for 
a particle outside a cylindrical shell as a function oi a/R for 
R/Xa = 1, 2.5 and 5 (solid lines 1, 2, and 3, respectively). 
The dashed line represents the Casimir-Polder potential in the 
static limit. 



1.1 




0.5 1.0 1.5 2.0 

a/R 



Fig. 3. Comparison between the exact and PFA Casimir- 
Polder potentials for particle-cylinder interaction (the solid and 
dashed lines, respectively) for R/Xa = 1 and 5 (lines labeled 
1 and 2) in the separation region 0.1 < a/R < 2. The inset 
shows the region from a/R — 0.01 to a/R — 0.1. 



different in nature. Particularly, the van der Waals force 
does not depend on c whereas the potential (|28l) in the 
limit a/R — > does. In fact, the nonrelativistic limit can- 
not be achieved for an ideal metal cylinder. Figure [5] shows 
the ratio [/(^*^')/{7 as a function of a/R (the upper, mid- 
dle and lower lines correspond to R/Xa — 1, 2.5, and 5, 
respectively) . 

Now we turn to the comparison between the exact 
Casimir-Polder energy with that obtained using PFA. The 



latter is given by the expression 



f/pFA(a) 



he 



R 



2C + C - 



R + a 
a(l4 



2(i?- 



(31) 



obtained from corresponding equation in |25| by putting 
the reflection coefficients equal to ±1, as it holds for ideal 
metal. In the static limit, (|3T|) leads to 



2.2 - 



2.0 




0.5 1.0 1.5 2.0 

a/R 

Fig. 2. The ratio of exact Casimir-Polder potentials calcu- 
lated using static and dynamic polarizabilities of a particle as 
a function of a/R for the same values of R/Xa, as in Fig. [U 
increasing from top to bottom lines. 



and for a/R ^ 1 one has 

t/pFA(a) « Uoia) (1 - 1^) - (33) 

to be compared with the exact result (|30p obtained in this 
limiting case. In Fig. [3] we plot the ratio U/Uq evaluated 
in the single-oscillator model by using exact expression 
(solid lines) and by the PFA, expression ([31]) (dashed 
lines). The lines 1 and 2 correspond to R/Xa — 1 and 5, 
respectively. As can be seen in Fig. [31 the PFA results are 
in very good agreement with the exact results. For exam- 
ple, at a/R = 0.1 the relative deviation between the two 
results is equal to 1.3% and 0.9% for lines 1 and 2, respec- 
tively. The same lines within the region 0.01 ^ a/R ^ 0.1 
are shown on an inset to Fig. [3] At a/R = 0.02 the relative 
deviation between the exact and the PFA results does not 
exceed 0.4%. 
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6 Microparticle in an interior of cylindrical 
shell 

Here, we present a few computational results for a mi- 
croparticle inside an ideal metal cylindrical shell. Substi- 
tuting (27) into (25) and integrating over ^, we find the 
expression 



Uir) 



he 



E 



FmiM 



df3 p 



,(/3)[s,(/3) + l] /„(/3i?) 



ifniM , F„ 



(34) 



(/3r) 



In the static limit, similar to (29), one obtains 



Krnim 
" Irnim) 



[F,„(/3r)+2/2,(/3r)] 



(35) 



The computational results are presented in Fig. 21 where 
the exact Casimir-Polder energy p4l) normalized for Uq is 
plotted as a function of a/i? (note that in the interior of a 
cylindrical shell a = R — r holds). In so doing the single- 
oscillator model was again used. Similar to Fig. [U the 
solid lines from the bottom to the top one correspond to 
R/^A = 1, 2.5, and 5, respectively. The dashed line repre- 
sents the Casimir-Polder potential in the static limit (I55|) . 
As is seen in Fig. |4j with increasing R/Xa the Casimir- 
Polder energy for a particle inside a cylindrical shell ap- 
proaches the static limit. This is in analogy to a particle 
outside an ideal metal cylindrical shell. 



7 Conclusions and Discussion 

In the foregoing we have derived the exact Casimir-Polder 
potential for a polarizable particle placed in an exterior or 
in an interior of an ideal metal cylindrical shell. Derivation 
was performed by the Green tensor method. The obtained 
results for a particle outside the cylindrical shell were com- 
pared with obtained by another method and found to be 
in agreement. Computations were made for particles de- 
scribed by the single-oscillator model. It was shown that 
at particle-cylinder separations below O.IR errors intro- 
duced by the use of the PFA do not exceed 1.3%. This 
justifies the application of the PFA for the comparison of 
the experimental data with theory. We have also inves- 
tigated the role of dynamic polarizability of the particle. 
According to our computational results, the static limit is 
approached for large a/R or R/Xa- Replacing ideal metal 
boundary conditions used in expressions (|6]) and Q by, 
for instance, standard electrodynamic continuity condi- 
tions or boundary conditions of the Dirac model pSl[?7] , 
it is possible to find the respective exact solutions for a 
particle interacting with a dielectric cylindrical shell or 
with a single-walled carbon nanotube. 

The obtained results can be used in the theory of topo- 
logical defects, e.g., when there is a cosmic string along 
the axis of the cylinder. For a cosmic string with the an- 
gle deficit 2tt — (po (0 ^ cj) ^ (po), the formulae for the 
Casimir-Polder potential are obtained from (fTO)) and 
by the replacement m — qui and adding the factor q, 
where q = 27r/0o (for the corresponding eigenfunctions 
see [2S^). Note that the conducting cylindrical surface can 
be considered as a simple model of superconducting string 
core. Superconducting strings are predicted in a wide class 
of field theories and they are sources of a number of in- 
teresting astrophysical effects such as generation of syn- 
chrotron radiation, cosmic rays, and relativistic jets. 
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Fig. 4. The normalized exact Casimir-Polder potential for 
a particle inside a cylindrical shell as a function of a/R for 
R/Xa = 1, 2.5 and 5 increasing from bottom to top solid lines. 
The dashed line represents the Casimir-Polder potential in the 
static limit. 
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